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Abstract 



o 

QJ ' The observation of duality by Mukherji and Mishra in one dimensional transport problems has been used to 

develop a general approach to classify and characterize the steady state phase diagrams. The phase diagrams are 
determined by the zeros of a set of coarse-grained functions without the need of detailed knowledge of microscopic 

jrt I dynamics. In the process, a new class of nonequilibrium multicritical points has been identified. 

-)— > ' 

^ ; 1 Introduction 

■ . ' There are many situations that involve transport of particles from one end to other along a one dimensional track 

'^ , obeying some form of mutual exclusion^]. Examples are vehicular traffic in a one-lane road, molecular motors 
fi ■ carrying cargo on a track in biological systems and so on. Moreover, such simple systems arc of importance to 
develop an understanding of nonequilibrium steady state phases and phase transitions. 

Many such transport models have now been studied extensively via analytical methods, meanfield approximations 

and numerical simulations, and different types of phases have been identified|21|31 01[S1|H1 |S]. These phases are 

represented in phase diagrams in the space of the externally controllable parameters of the problem, namely, the 

imposed rates of injection and withdrawal at the two boundaries. The fact that the boundary conditions determine 

^+ ' the stable phase diagrams makes these nonequilibrium problems different from equilibrium ones. 

^^ , Since phase transitions involve the whole of the system, the generic behavior, in the large size limit, is expected 

\l ■ to be determined by certain gross overall features. This is the way equilibrium phase transitions are analyzed, but, 

alas, no such general framework is known for nonequilibrium cases. Hence the efforts on case by case studies. Our 

aim is to develop a general formulation at least for the above class of systems. We show that the generic features 

and the universal properties of the phase diagrams can be obtained from the structure of a set of S'-functions and 

C^ ' by using a recently discovered duality |3j, without the need of detailed specifics of the microscopic dynamic rules and 

G , interaction. Still, one microscopic parameter remains essential for the problem, namely, a small distance cutoff (e.g. 

lattice spacing or some microscopic size etc) which cannot blindly be set to zero. The usefulness of the approach is 

shown by predicting a new class of nonequilibrium multicritical points. 



2 Phases and response functions 



k>( I As an example 12 El El J consider a one dimensional lattice of iV(— >• oo) sites. Particles are injected at site i = 1 at 
; I ' a rate a (i.e probability that a particle is injected in a short time interval dt is adt) and withdrawn at site i = N 
Cd I at a rate 1 — 7. The particles hop on the lattice as per preassigned rules, like forbidden multiple occupancy of a 
site, etc. In addition, non-conserving processes may allow addition to or deletion from the track at rates Ua and ujd 
respectively, as e.g., exits or feeders in a traffic system for cars to get out of or into the road, or "processive" objects 
in biological systems falling off the track or getting reattached from the bulk solution. The parameters, K = uJa/i^d 
and n = ujdN are characteristics of the microscopic dynamics while a and 7 are externally imposed. 

For a coarse-grained description, the natural variable is the local density or the space-time dependent average 
occupation number p{x, t) in continuum {x G [0, 1] by a rescaling of the total length). The sensitivity to the boundary 
concentrations (or rates) can be measured by the response functions 

Xfj. = -t: — , where /i = a or 7, and AI = dx p{x), (1) 

C'M Jo 



is the steady state spatially averaged density. Any two points in the — 7 space are said to be in the same phase 
if they can be connected by a continuous path along which the density profile or the response functions change 
smoothly. Any point of non-analyticity on a path defines the location of the phase transition. 

The plrascs observed are, generally, of the following types, (i) Injection (withdrawal) rate dominated, to be called 
the a-phase (7-phase), (ii) a shock phase consisting of piecewise continuous densities, and (iii) special phases. In the 
shock phase, there is a discontinuity in the density separating an a-phase on one side from a 7-phase on the other 
side, while an example of (iii) is a phase where tire current through the system is maximum. The response functions 
behave differently in these phases. In the a-phase, Xa ~ 0(1) but x-y = Oi while in the 7-phase, it is the other way 
round. However, in the shock phase, both Xa and X7 would be nonzero. For special phases like (iv) in the above 

list, Xa,7 ~ 0. 

3 Equations for dynamics and steady states: Definitions of Si's 

In a continuum limit for large N (with lattice spacing a -^ L = No) the time variation of p[x, t) can be written in 

the form of a continuity equation, as 

dp{x) dJoix,t) 

^^ + ^^=5o(p,i), (2) 

where the right hand side is the explicit non-conserving contribution to the change in density. The left hand side 
is in the form of a continuity equation with Jo(a;,t) as the current at the site. In a mean field approximation, the 
current is taken to be an implicit function of position and time through the density, so that Jq can be split into two 
parts, 

Jo = -eS2{p)^ + J{pix,t)), (3) 

with a "bulk" contribution J(p) determined by the local density and a term that depends on the derivative of the 
density ('Tick's law") over the lattice spacing, e = a/2N being small. This e-dependent term is a reminiscent of 
the interactions in the neighborhood of a site on the underlying lattice. The form of 5*2 (p) is determined by the 
microscopic dynamics, but, for simplicity, we take 52 (p) = 1 here. Such a form like Eq. |31has recently been shown 
by Chakrabarti[Sj to arise naturally in a renormalization group type approach for transport processes and failures of 
fiber bundles. The fact that there are two fixed points (viewed as a recursion relation) will be of importance to us 
also. 

In the steady state, the system evolves to a time independent density profile satisfying 

-e^S2{p)^ + S,{p)^ + So{p)^0. (4) 

ax ax ax 

The Si functions encode the dynamics or specialties of the system. The density satisfies the boundary conditions 
p{0) = a and p{l) =7. The microscopic rules are taken to be sufficiently smooth to warrant considerations of smooth 
and analytic functions only. These restrictions can be relaxed if necessary. 

For e — » 0, the ensuing first order equation cannot in general satisfy the two boundary conditions. Therefore, the 
e term, eventhough looks innocuous in the bulk limit, is essential. It defines a new scale x = x/e in the problem and 
this scale is important for the phase transitions. E.g., the discontinuity at a shock will be rounded on a scale of x 
but would look sharp on a bigger scale. 

3.1 Nature of >So 

Let us first consider the role and the nature of So{p). The zero of the non-conservation function So{p) is a special 
density. This is the equilibrium like steady state density p ~ p^, Sq{pl) = 0, the system would evolve to if all other 
dynamics, except this non-conserving one, are switched off. In such a situation, the density pL can be obtained from 
the extrema of a (free energy like) Lyapunov function, So{p) such that So{p) = —dSo/dp. For stability of the state, 
evaporation is to be preferred in case of excess density (over p^), but adsorption for p < p^. This requires Sq{p) to 
be an odd function or, Sq{p) an even function oi p ~ pL- A simple possible choice is 

So{p)^nip-pLf', iq>l), (5) 

with q = 1 corresponding to a linear form for 5*0 (the so-called Langmuir kinetics). The "softness" of the state is 
determined by the parameter il that controls the width of the well. Also, the conserved case is recovered by taking 
r2 — > 0. A bistable (or multi-stable) situation can be obtained for fi < with additional terms in So. 



3.2 Nature of Si 

A zero of Si{p), i.e. Si{pm) = 0, is the density at which the current is an extrema (e.g. a maximum). From 
Eq. Q, we see that this is the density where one may afford a non-existence of the first derivative of the density. 
Consequently, a shock, for which the derivative is not defined (strictly in the e — > limit), if exists, has to be centered 
around p ~ p„j. If the dynamics has a particle hole symmetry, then p„i = 1/2. For concrctencss and simplicity we 
consider the class of functions 

J(p) « J„, - a(p - p„0'^ (P>1), (6) 

near the maximum. There can be cases with more than one zero of Si{p), which can lead to multiple shocks and 
more exotic phenomena. Such cases will be discussed elsewhere. 

3.3 Special cases 

The microscopic dynamic rules determine the values of the two special densities, pL,Pm, and the values oi p,q. 
However, we do not require those rules. We need to distinguish special cases like, (i) pL > Pm, (h) PL < Pm,and (iii) 
Pl = Pm- It is plausible to make a smoothness hypothesis that the nature of the phase diagram in the a~"f (external)- 
parameter space changes smoothly as the parameters like pi , p„j (determined by the microscopies) are changed unless 
there is a special condition. Such a condition is p^ — Pm where the non-conserving processes try to maintain a 
density at which the conserved processes can accommodate maximum current. 

4 Boundary layer approach: outer and inner densities 

We adopt the boundary layer approach or the method of matched asymptotics to handle the two scales in Eq. 10}. 
Consider the case where the bulk density profile p — Pontix), obtained from Eq. ^with e = 0, matches the boundary 
condition p{x = 0) = a. But, then, po = Pont{x = 1) 7^ 7. A different density profile Pin{x) where x = {x — Xs)/e 
^ 0(1), {xs = 1), extrapolates within a thin "inner" region from (Oin(— 00) ~ po to Pin{0) = 7- This inner region 

satisfies (see Eq^ Si{p) = dSi{p)/dp] 

-52(pi„)^ + 5i(pi„)=0, (7) 

ax 

which is equivalent to Eq. |21with Si{p) = J(p) — Jq. The inner region, to first order in e, is too thin for the violation 
of conservation to matter, so that the current Jo = J{po) entering from the bulk (outer) region remains conserved 
in the inner layer. Inference: The mandatory matching condition requires Si{p) to have a zero at p ~ po- 

4.1 Zeros as requirements for shocks 

A shock is formed only if the inner solution fails to satisfy the boundary condition. This happens if the inner solution 
saturates at the other end. Therefore the minimal requirement for shock formation is another zero of Si (p) , so that 

Si{p) = -{p-po){p-PsMp). (8) 

The first nontrivial case, is therefore a function with two simple zeros and $(p) — 1. The two zeros Po and Ps 
correspond to the two fixed points of Chakrabarti's approach|5]. By RoUe's theorem, po < Pm < Ps- 

The inner equation admits two types of solutions, one bounded (B-type) between po and ps while the other one 
shows a divergence (U-type) with dp/dx ^ —p^, or more generally, dp/dx ~ —p^^- It is the B-type layers that 
matures to a shock but not the U-type. The inner solution is of the form T(x/w + C) with I{z) — > ps or po as 
z -^ ±00. Here w is the width of the layer and ^ gives the location of the center of the layer. So instead of the 
two boundary conditions describing the layer, we may instead opt for the w,^ pair. The center may lie outside the 
physical range or may be in an unphysical density range, requiring continuation of the density and the space beyond 
the physical range of [0, 1]. This continuation helps in getting the general form of the phase diagram. The origin is 
to be called a "virtual origin" if it is in the unphysical region. 

4.2 Shockening transitoin and Mukherji-Mishra dual line 

For a given a, as 7 is changed, two different situations may arise. In one, the virtual origin approaches the boundary 
at X = (i.e. x = 1) and then enters the physical region, eventually moving to —00. In the other situation, the 
origin remains virtual and moves to infinity, ^ — *■ -l-oo. This is the Mukherji-Mishra (MM)dual boundary line. 



The first case is a tliickcning of tlie layer but remaining pinned to the boundary. Ultimately as ^ -^ — cxo, the 
layer gets released from the boundary and moves into the bulk. Or a shock forms. So long as the boundary layer 
stays pinned to the boundary, x-y ~ ^l/^i{l) —^ as e ~* 0. In contrast, Xa is nonzero. The phase, by definition, is 
then an a-phase. 

The transition of a thin layer to a shock at 7 = ps{po{ct)) has been called a "shockening" transition or a layer 
"shockens". Beyond this, the layer is separated from the boundary by a bulk phase (outer solution) of nonzero 
thickness. Though the shockening of the inner layer is a depinning phenomenon at the boundary, it also signals a 
bulk phase transition from an a-phase to a shock phase. It is apparent that the shock phase has both the response 
functions Xq,7 ¥" 0- 

The symmetry of the two zeros of Si{p) suggests that there has to be another line 7 = Poict) at which ^ —> +00. 
The boundary region goes from an accumulated to a depleted region as one crosses this MM line, thereby separating 
the shockening (B-type) to nonshockening (U-) type boundary layers. The MM dual line is purely a boundary 
transition line, and its existence is a requirement for shock formation. 



4.3 Two lengthscales 

For 7 near the two extreme values poj Psj the lengthscale ^ can be obtained as ^ ~ X~^('y/ px) where X stands for 
o or s as appropriate and /^^ is the inverse function of X (defined for the inner solution). Eq. IJHJ) suggests, T~^ 
to be logarithmic implying ^ ~ log \ j ~ px I • We note here that from the exact solution of the totally asymmetric 
exclusion problem (with conservation), one can associate this dual line {a ~ 7) with ^ ^ log | a — 7 |, identical to 
the result we just derived. 

(b) 

7 




MM dual line 








\ '^ ' 




■«■ ■ 


\ 


/ / 


■^x 


■H« .- 




,-■ 










%v 






Figure 1: Possible phase boundaries (red thick line) and dual lines(blue dotted). (a,b) for the a-phase with the shock 
forming at .t = 1 while (c,d) for the 7-phase with shock forming at a; = 0. In (a) and (c) the intersection of the shock 
phase boundary and the MM dual line produces a critical point marked by a filled square. No critical point in (b,d). 
The a, 7 asymmetry depends on the relative values of pL and pc- 

The other length scale w can be obtained from various limits of Eq. Q, the lengths differing by constant 
factors. From the asymptotic approach to the limits px, (X=o or s), w^^ = {ps — Po)^{px), while for p w p,,, one 
gets w~^ ^1 Po — Ps I /Si[pm)- What is important to note is that for a given a, the width is determined by the 
corresponding separation of the shockening and the dual line, to be called the MM-gap. The height h of the shock 
on the shockening line is also equal to this MM-gap. 

4.4 Condition for Critical point 

In case the shockening transition line and the dual line intersect, then the intersection is at 7 = pc with iti — !• 00 as 

w^\p,^p,\-'^^P~^^^h-'^^P-^\ (9) 

Such a divergence is the signature of a critical point. The bulk phase transition from the a-phase to the shock phase 
is first order because at the transition point h > 0. On the other hand, the shock evolves from a zero height at the 
critical point so that it is a continuous transition. In case the two lines do not cross, there will be no critical point 
and the lines will span the whole phase diagram, symmetrically placed around 7 = pm if *I'(p) = 1- 



4.5 Phase diagrams 

So far we have concentrated on the a-phase only. A similar analysis can be done for the 7-phase for which the 
shock is formed at x = 0. Here again there are two possibilities; the shockening and the dual lines either intersect 
at a = Pm or do not intersect but remain on two sides of a = pm- AH the four possibilities are shown in Fig. ^ In 
these diagrams the lines at a = ac or 7 = 7c or both remain special like the dual lines, representing boundary layer 
transitions. 
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Figure 2: a — 7 phase diagrams. (a,b,c) for p^ = Pm, (c) for p^ > pm and (d) for pL < pm- The critical points 
are represented by filled circles. The shaded region is the shock region. In (a) q > p. (b) q = p while (c) has no 
non-conserving part. The MM dual lines are shown in (c) and (d) only. The filled triangles on the dual lines represent 
the point {pl,Pl) which is always on the dual line. 

Combining the two, we can now draw the global phase diagram. Combination of (a) and (c) of Fig. ^one gets 
the type known for the pi = p,n [K = 1) case of Ref. |21E1IH1- Similarly, (a) with (d) is known for p^ > Pm, 
while (b) with (c) will be the case for p^ < Pm- For 7 < 7c = Pm, on the a-phase side, the U-type boundary layer 
renders an effective boundary value pm at a; == 1 and therefore the critical behavior continues for all 7. The shock on 
changing a evolves in height and shifts to x = 0. On the 7-phase side, the response function x-y undergoes a change 
on crossing the line 7 = pm, even though the bulk density distribution changes smoothly. The line 7 = pm, like the 
MM dual lines, indicates a boundary transition, but in special situations it may develop into a bulk phase boundary 
also. The latter happens as p^ — > Pm- 

The shape of the shockening curve near the critical point for (7 = 1 is given by (7 — 7^) '^l a ~ ac 1^"^. For 
7 = 7c = Pm, the shock height vanishes on the shock side a.s h ^\ a — ac |^ with C,' = l/{2p+ 1). Though the 
p = 1 case is known in the literature, we have identified the whole class of multicritical points. All the exponents^] 
associated with the shockening transition and the critical point can be determined in terms of g, p. These details will 
be reported elsewhere. 

For PL = Pm, our analysis via the duality yields the nature of the critical points. For q = p, the phase boundaries 
are similar to the q = p = 1 case. However for q > p, the critical point is at a = p„,, 7 = pm- We show a new type 
of phase diagram for a particular case with pm ~ 0.5 in Fig. |21 

As one traverses the shock phase from one phase boundary to the other, the shock position goes from x = 1 to 
X = 0. Now, if the non-conserving part of the dynamics is removed, the shock region collapses on to a line as in 
Fig. I21c). But the collapse also means that the shock is uniformly distributed over the entire length and the density 
is to be averaged over this distribution of shocks [TO]. This yields the linear density profile one knows from exact 
solutions^]. This also shows that the mean field theory puts a bias towards shock formations, so that a judicious 
use is called for in situations where shocks are not expected. 

5 Summary 

In summary, the MM duality theorem can be stated as follows: (a) Every shockening transition has a dual boundary 
transition, (b) If the two lines (the shockening transition and the dual line) intersect, there is a critical point, (c) 
The nature of the critical point is determined by the zeros of the S'^-functions. This theorem has been used to predict 
the behavior of a class of multicritical points in the steady state phase diagram of nonequilibrium transport. Our 
conclusion is that all microscopic perturbations need not be relevant to the phase diagrams. The microscopic dynamics 
needs to be analyzed for the nature of the zeros of the Si functions and the values of g, p. Those interactions or rules 
that change the density parameters pL,m without change in q^p, can be grouped into the same class. Perturbtions 
within the class will only make cosmetic changes in the phase diagram. We have shown a few examples of possible 
multicritical phase diagrams. New classes can be generated by including extra features of 5'i's. 
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